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Global analysis of neutrino oscillation data in four-neutrino schemes 
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We present an analysis of the global neutrino oscillation data in terms of four-neutrino mass schemes. We find 
that the strong preference of oscillations into active neutrinos implied by the latest solar as well as atmospheric 
neutrino data allows to rule out (2-f 2) mass schemes, whereas (3-1-1) schemes are strongly disfavoured by short- 
baseline experiments. Our analysis shows that four-neutrino oscillations do not provide a satisfactory description 
of the global neutrino oscillation data including the LSND result. 



1. INTRODUCTION 

The neutrino oscillation interpretations of the 
solar and atmospheric [||J^ neutrino data 

and the LSND experiment ^ require three neu- 
trino mass-squared differences of different orders 
of magnitude. Since it is not possible to ob- 
tain this within the Standard Model framework 
of three active neutrinos it has been proposed 
to introduce a light sterile neutrino to rec- 
oncile all the experimental hints for neutrino os- 
cillations. Here we present an analysis of the 
global neutrino oscillation data in terms of four- 
neutrino mass schemes, including data from solar 
and atmospheric neutrino experiments, the LSND 
experiment, as well as data from short-baseline 
(SBL) experiments |]7|,p|,p| and long-basehne re- 
actor experiments pOl reporting no evidence for 
oscillations. This analysis updates the work pre- 
sented at the meeting and is based on the data of 
summer 2002. We find that for all possible types 
of four-neutrino schemes different sub-sets of the 
data are in serious disagreement and hence, four- 
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neutrino oscillations do not provide a satisfactory 
description of the global oscillation data includ- 
ing LSND. The details of our calculations can be 
found in Refs. 

Four-neutrino mass schemes are usually di- 
vided into the two classes (3-1-1) and (2-f2), as 
illustrated in Fig. |l|. We note that (34-1) mass 
spectra include the three-active neutrino scenario 
as limiting case. In this case solar and atmo- 
spheric neutrino oscillations are explained by ac- 
tive neutrino oscillations, with mass-squared dif- 
ferences Am^QL ^-nd Ar7i^.j,^,, and the fourth neu- 
trino state gets completely decoupled. We will 
refer to such limiting scenario as (3-1-0). In con- 
trast, the (2-1-2) spectrum is intrinsically differ- 
ent, as there must be a significant contribution 
of the sterile neutrino cither in solar or in atmo- 
spheric neutrino oscillations or in both. 



2. NOTATIONS AND 
TIONS 



APPROXIMA- 



Neglecting CP violation, in general neutrino os- 
cillations in four-neutrino schemes are described 
by 9 parameters: 3 mass-squared differences and 
6 mixing angles in the unitary lepton mixing ma- 
trix. Here we use a parameterization introduced 
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Figure 1. The six types of four-neutrino mass 
spectra, divided into the classes (3+1) and (2+2). 



in Ref. which is based on physicaUy rele- 

vant quantities: the 6 parameters ArngQ^i ^sol, 
Am^^^,, 61 ATM, AttIlsnd' ^lsnd are similar to the 
two-neutrino mass-squared differences and mix- 
ing angles and are directly related to the oscil- 
lations in solar, atmospheric and the LSND ex- 
periments. For the remaining 3 parameters we 
use ris,rie and d^. Here, ri^ [rje) is the fraction 
of Ps {ve) participating in solar oscillations, and 
(1 — d^) is the fraction of i^^ participating in os- 
cillations with Ato^^j, (for exact definitions see 
Ref. For the analysis we adopt the follow- 

ing approximations: 

• We make use of the hierarchy 

« ^rnlru « ^mlsNu ■ (1) 

This means that for each data set we con- 
sider only one mass-squared difference, the 
other two are set either to zero or to infinity. 

• In the analyses of solar and atmospheric 
data (but not for SBL data) we set rye = 1, 
which is justified because of strong con- 
straints from reactor experiments 

Due to these approximations the parameter 
structure of the four- neutrino analysis gets rather 
simple. The parameter dependence of the differ- 
ent data sets solar, atmospheric and SBL is illus- 
trated in Fig. H We see that only r]s links solar 
and atmospheric data and links atmospheric 
and SBL data. All the other parameters are "pri- 
vate" to one data set. 




Figure 2. Parameter dependence of the different 
data sets in our parameterization. 



3. (2+2): RULED OUT BY SOLAR AND 
ATMOSPHERIC DATA 

The strong preference of oscillations into ac- 
tive neutrinos in solar and atmospheric oscilla- 
tions leads to a direct conflict in (2+2) oscilla- 
tion schemes. We will now show that thanks to 
the new SNO solar neutrino data |^ and the im- 
proved SK statistic on atmospheric neutrinos |^ 
the tension in the data has become so strong that 
(2+2) oscillation schemes are essentially ruled 
outj^ Indeed, latest solar neutrino data lead to 
the bound rjs < 0.45 at 99% CL, where 77^ is the 
parameter describing the fraction of the sterile 
neutrino participating in solar neutrino oscilla- 
tions. In contrast, in (2+2) schemes atmospheric 
data imply rjs > 0.65 at 99% CL, in clear dis- 
agreement with the bound from solar data. 

In Fig. ^ we show the for solar data and for 
atmospheric combined with SBL data as a func- 
tion of rjs- Furthermore, we show the of the 
global data defined by 

X^iT^s) = AxL(^s) + Axl^,„iVs) ■ (2) 

From the figure we find that only if we take both 
data sets at the 99.95% CL a value of rjs exists, 
which is contained in the allowed regions of both 
sets. This follows from the x^-value Xfc ~ -'^2.2 
shown in the figure. In Ref. we have proposed 
a statistical method to evaluate the disagreement 
of different data sets in global analyses. The pa- 
rameter goodness of fit (PG) makes use of the x^ 

^Details of our analyses of the snl ar a .nri atmospheric 
neutrino data can be found in Ref. ]l3| . For an earlier 
four-neutrino analysis of solar and atmospheric data see 
Ref. 0. 
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Figure 3. AxgoL> ^xItm+sbl and Xgiobai as a func- 
tion of in (2+2) oscillation schemes. 



defined in Eq. (^. This criterion evaluates the 
GOF of the combination of data sets, without be- 
ing diluted by a large number of data points, as it 
happens for the usual GOF criterion (for details 
see Ref. @). We find xIg = xLn = 23.5, lead- 
ing to the marginal PG of 1.3 x 10^^. We con- 
clude that (2-f 2) oscillation schemes are highly 
disfavoured by the disagreement between the lat- 
est solar and atmospheric neutrino data. This is a 
very robust result, independent of whether LSND 
is confirmed or disproved. 

4. (3+1): STRONGLY DISFAVOURED 
BY SBL DATA 

It is known for a long time that (3+1) 
mass schemes are disfavoured by the comparison 
of SBL disappearance data with the LSND 
result. In Ref. ||l^ we have calculated an up- 
per bound on the LSND oscillation amplitude 
sin^ 20LSND resulting from SBL and atmospheric 
neutrino data. From Fig. ^ we see that this 
bound is incompatible with the signal observed 
in LSND at the 95% CL. Only marginal over- 
lap regions exist between the bound and LSND if 
both are taken at 99% CL. An analysis in terms of 
the parameter goodness of fit ||l^ shows that for 




Figure 4. Upper bound on sin^ 20lsnd from SBL 
and atmospheric neutrino data in (3+1) schemes 
compared to the LSND allowed region ||]. 



most values of Am^gp^^ SBL and atmospheric data 
are compatible with LSND only at more than 
3(7, with one exception around Amlg^^ ~ 6 eV^, 
where the PG reaches 1%. These results show 
that (3+1) schemes are strongly disfavoured by 
SBL disappearance data. 

5. COMPARING (3+1), (2+2) AND 
(3+0) HYPOTHESES 

With the methods developed in Ref. ||T^ we 
are able to perform a global fit to the oscilla- 
tion data in the four-neutrino framework. This 
approach allows to statistically compare the dif- 
ferent hypotheses. Let us first evaluate the GOF 
of (3+1) and (2+2) spectra with the help of the 
PG method described in Ref. ||l^. We divide the 
global oscillation data in the 4 data sets SOL, 
ATM, LSND and NEV, where NEV contains the 
experiments KARMEN §, GDHS §, Bugcy § 
and CHOOZ/Palo Verde ||l^, reporting no evi- 
dence for oscillations. Following Ref. ||l^ we con- 
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SOL ATM LSND NEV 


XPG 


(3+1) 
(2+2) 


0.0 0.4 7.2 7.0 
14.8 6.7 2.2 9.7 


14.6 
32.4 



Table 1 



Contributions of different data sets to Xpc ™ 
(3+1) and (2+2) schemes. 



sider 

= Axsol(^sol, Am2oL>'7s) 

+ AXatm(^atm, Amj^.j,^,[, 77s, d^) , , 

+ AXnev(^lsnd, ATO2s„o,d^,7?e) 

where ^ x\ ~ (x^)min {X = SOL, ATM, 

NEV, LSND). In Tab. g we show the contribu- 
tions of the 4 data sets to = Xmin fo^' (3+1) 
and (2+2) oscillation schemes. As expected we 
observe that in (3+1) schemes the main contri- 
bution comes from SBL data due to the tension 
between LSND and NEV data in these schemes. 
For (2+2) oscillation schemes a large part of Xpg 
comes from solar and atmospheric data, however, 
also SBL data contributes significantly. This 
comes mainly from the tension between LSND 
and KARMEN 0, which does not depend on 
the mass scheme and, hence, also contributes in 
the case of (2+2). Therefore, the values of Xpq i^i 
Tab. |l| for (2+2) schemes are higher than the one 
given in Sec. |3[ where the tension in SBL data is 
not included. 

The parameter goodness of fit is now obtained 
by evaluating Xpc f^i' 4 DOF jl2j. This number 
of degrees of freedom corresponds to the 4 pa- 
rameters rjsjdfj^, 6'lsnd, AttIlsnd describing the cou- 
pling of the different data sets (see Eq. (||)). The 
best GOF is obtained in the (3+1) case. How- 
ever, even in this best case the PG is only 0.56%. 
The PG of 1.6 X 10"^ for (2+2) schemes shows 
that these mass schemes are essentially ruled out 
by the disagreement between the individual data 
sets. 

Although we have seen that none of the four- 
neutrino mass schemes can provide a reason- 
able good fit to the global oscillation data in- 
cluding LSND, it might be interesting to con- 
sider the relative status of the three hypothe- 
ses (3+1), (2+2) and the three-active neutrino 



scenario (3+0). This can be done by compar- 
ing the x^ values of the best fit point (which is 
in the (3+1) scheme) with the one correspond- 
ing to (2+2) and (3+0). First we observe that 
(2+2) schemes are strongly disfavoured with re- 
spect to (3+1) with a Ax^ = 17.8. The reason 
for the big change with respect to the value of 
Ax^ =3.7 found in Ref. jl|] is the improved sen- 
sitivity of solar (SNO NC) and atmospheric (SK 
1489-days) data to a sterile component. With this 
new data now (3+1) schemes are clearly preferred 
over (2+2): for 4 DOF a Ax^ = 17.8 implies that 
(2+2) is ruled out at 99.87% CL with respect to 
(3+1). Further we find that (2+2) is only slightly 
better than (3+0), which is disfavoured with a 
Ax^ = 20.0 with respect to (3+1). 

6. CONCLUSIONS 

Using latest solar and atmospheric data, we 
re-analyze the four-neutrino description of cur- 
rent global neutrino oscillation data, including 
the LSND evidence for oscillations. The higher 
degree of rejection for non-active solar and atmo- 
spheric oscillation solutions implied by the SNO 
neutral current result as well as by the latest 
1489-day SK atmospheric neutrino data allows us 
to rule out (2+2) oscillation schemes. Using an 
improved goodness of fit method especially sensi- 
tive to the combination of data sets we obtain a 
GOF of only 1.6 x 10"^ for (2+2) schemes. Fur- 
ther, we find that also (3+1) schemes are strongly 
disfavoured by the data. The disagreement be- 
tween negative short-baseline experiments and 
LSND in (3+1) schemes imply a GOF of 5.6 x 
10"'^. This leads to the conclusion that all four- 
neutrino descriptions of the LSND anomaly, both 
in (2+2) as well as (3+1) realizations, are highly 
disfavoured. Our analysis brings the LSND hint 
to a more puzzling status, and the situation will 
become even more puzzling if LSND should be 
confirmed by the up-coming MiniBooNE experi- 
ment 
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